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1 Introduction 



Several measurements have reported evidence for flavor mixing in the neutral D me- 
son system [TJ EJ [3j HJ El E] . The results are in agreement with the Standard Model 
(SM) predictions, which unfortunately are affected by large theoretical uncertain- 
ties [7J EJ [9] . However, the increasing precision of D°-D° mixing measurements helps 
constrain new physics models pJJ [TT] . CP violation in D meson decays, though no- 
toriously difficult to calculate precisely, is expected to be very small in the SM, at 
the level of 10~ 3 or less [HI 121 EES]- Relatively large CP asymmetries, at the percent 
level, might be a signature of new physics effects. Recent results from the LHCb 
experiment [H] reported evidence for direct CP violation measuring the difference of 
CP asymmetries in singly- Cabibbo suppressed D — > tt + it~ - and D° K + K~ de- 
cays, with a statistical significance of 3.5a. The observed asymmetries are marginally 
compatible with the SM but not conclusive for establishing new physics [T5l [T6l ITT] . 
These intriguing results renew interest in studying mixing and CP violation in the 
D°D° meson system and in general in the charm physics sector. 

The Heavy Flavor Averaging Group (HFAG) [18] provides averages for heavy 
flavor quantities. It is divided into several sub-groups, each of which focuses on a dif- 
ferent set of heavy flavor measurements. The Charm Physics sub-group [JJJ] studies 
the following topics: D°-D° mixing, CP violation, spectroscopy of charm mesons and 
baryons, semileptonic decays, decay constants, hadronic branching fraction measure- 
ments and rare D decay modes. In particular, it provides averages for the the D°-D° 
mixing and CPV parameters by combining measurements from different experiments 
in a x 2_ based fit. The results of the fit are expressed in term of physics parameters 
that can be directly compared to the theoretical predictions. 

2 HFAG notations for flavor mixing and CPV 

Flavor mixing occurs when the Hamiltonian eigenstates (or physics eigenstates) Di, 
D 2 differ from the flavor eigenstates D°, D°. A neutral D meson produced at time 
t = in a definite flavor state D°, will then evolve and oscillate into a state of 
opposite flavor D°, after a certain time. If we describe the time evolution of the 
flavor eigenstates in terms of a 2 x 2 effective Hamiltonian, H = M — |T, then, 
assuming CPT is conserved, the mass eigenstates can be expressed in terms of the 
flavor eigenstates by, 

| J D 1 )=p| J D ) + g|D°) (1) 
\D 2 ) = p\D°) - q\D } 

*The use of charge conjugate reactions is implied throughout unless stated otherwise. 
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with the normalization \q\ 2 + \p\ 2 = 1 and 
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Assuming a phase convention such that CP\D°) = —\D°) and CP\D°) = —\D°) then, 
if CP is conserved, we have that q = p = l/\/2 and the mass eigenstates coincide 
with the CP eigenstates: |P>i) = |P>cp_> (CP-odd) and |D 2 ) = |P>cp + > (CP-even). 

The mixing parameters can be expressed in terms of the difference of masses (7^1,2) 
and widths (I\ 2 ) of the Hamiltonian eigenstates, 



x 



TX12 — m\ 



2r 



(3) 



where T = (T 1 + T 2 )/2. 

CP violation can be of three types: 



CPV in decay or direct CPV: this occurs when the decay amplitudes for CP 
conjugate processes are different in modulus. If (f\H\D°) — Af, (f\H\D°) = Aj 
are the D° and D° decay amplitudes into the final states / and CP conjugate 
/, then 

A r 

± 1 =^ CPV. (4) 



2. CPV in mixing or indirect CPV: it occurs when the Hamiltonian eigenstates do 
not coincide with the CP eigenstates. That is 



7^1 



CPV. 



(5) 



CPV in the interference of mixing and decay: for neutral D mesons there is a 
third possibility to observe CP violation even when CP is conserved in mixing 
and also in decay. In this case, CP violation arises when, in a process with final 
state / that can be reached by neutral D mesons of both flavors {i.e. D° and 
D°), there is a relative weak phase HI difference between the mixing and the decay 
amplitudes. The quantity of interest that is independent of phase conventions, 
and physically meaningful, is 



A 



/ 



q_Aj_ 

P A f 



A 



A 



(6) 



'The CP- violating phase is also indicated as a weak phase since it originates from the weak 
interaction in the SM. 
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where 5f and 0/ are the CP-conserving and CP- violating phases respectively. 
If CP is conserved in mixing and in decay, the signature of CP violation in the 
interference of mixing and decay is thus 



fa CPV. 

For CP eigenstates, CPV in either mixing or decay is indicated by 

|A/|^1, 

while CPV in the interference of mixing and decay corresponds to 

lm(X f ) ? 1. 



(7) 



(8) 



(9) 



Note that if there is no weak phase in the decay amplitudes then arg (q/p) = 4> 
and it is independent of the final state /. 



3 Experimental observables and parameters of the 
theory 

The most precise constraints on the mixing parameters x, y and the CPV parameters 
\q/p\, <f> are obtained in time-dependent analyses of D decays. 

Consider a final state / that can be reached by both D° and D° decays. If a D° 
is produced at time t — 0, it can reach the final state / by mixing to D° followed by 
the decay D° — > /, or directly through the decay D° — > f. The interference between 
the mixing and decay amplitudes modifies the time-dependence with respect to the 
pure exponential as follows: 



\(f\H\D°m 

Ip-rt 



lp-n 



= \A f \ 2 + 



+ 



+ 



\A f \ 2 + 



\A f \ 2 



\A f ? 



cosh(yrt) + 2Re \^A}A f j sinh(yn) 
cos(xTt) - 21m (^A* f A^J sm(xTt) (10) 



\A f \ 2 



cos 



h(yTt) + 2Re (^ A f A *^j sinh(yn) 



cos(xTt) - 21m I^A f A* f ) sm(xFt). (11) 



HFAG combines 38 observables measured in time-dependent and time-integrated anal- 
yses from the following experiments: BABAR, Belle, CDF, CLEO, CLEOc, E791, FO- 
CUS, and LHCb. When allowing for CPV there are 10 underlying parameters that are 
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extracted from the x 2 fit: x, y, \q/p\, 4>, 5, Sk™, Rd, Ad, A n , Ak- The parameters 
5, Sktttt are relative strong phases, R D is the ratio T(D° — y K + ir~)/T(D° — y K~ir + ) 
and Ad, A n , Ak are the direct CP-violation asymmetries for the D° — y K + ir~, 
D° — y 7r + 7r~ and D° —y K+K~ modes. The relationships between these parameters 
and the measured observables are given below. 

1. Semileptonic decays: search for mixing by reconstructing the "wrong-sign" (WS) 



decay chain, D*+ ->■ D°n + , P> u ->■ P> u , P> u ->■ 



#(*)- 



e zv In contrast to hadronic 



decays, the WS charge combinations can occur only through mixing. The mea- 
surement of Rm is related to the mixing parameters as follows 



Rm = -(x 2 + y 2 ), 



(12) 



and can be obtained directly as the ratio of WS to right-sign (RS) signal events. 
The RS events correspond to the non-mixed process. 

2. Decays to CP eigenstates: measure the mixing parameter yep and the CPV pa- 
rameter Ay with a lifetime ratio analysis of the transitions to the CP eigenstates 
and the transitions to the CP-mixed state D° — y K~n + , 
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The parameter A r is the decay-rate asymmetry for the CP eigenstates. If CP 
is conserved yep = y and A r = 0. 

Time- integrated CP- violating asymmetries in CP-even eigenstates, e.g. in D° — y 
K+K^ and D° — y tt+tt^ , provide constraints on the mixing and CPV parameters 
according to the relations, 
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(£) is the average reconstructed D° proper time and r D is the nominal D° 
lifetime. 
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3. Three-body D° — > fT°7r + 7r~ and D° — >■ K$K + K~ decays: measure directly the 
mixing and CiV parameters x, y, \q/p\, and with a time- dependent Dalitz 
plot analysis. 

4. Wrong-sign decays to hadronic non-CP eigenstates: measure the parameters 
x /=t and y f± and Re> and Ad in a time-dependent analysis of the WS events 
selected through the decay chain D* + — > D°% + , D° — > K + n~ . The parameters 



are defined as follows: 



1 ± A M \ 1/4 

(x' cos ± y' sin 0) (16) 



1/ = z -7— {y cos T a: sin 0) 

V J- =F %/ 

r(D° ->■ K+7T-) + T(D° ->■ K-7T+) _ 

r(£>° ir-7r+) + r(I7 K+7T-) ~ D 
r(D° ->■ x + 7r-) - r(_D° ->■ K-7T+) _ 
r(D° ^-tt+) + r(I7 fs:+7r-) ~ D 



where 





/ cos 5 sin S 
\ — sin 5 cos 5 

^ M k/p| 2 + b/?l 2 ' 6 arg 

The superscript in x f± and y /=t identifies the flavor of the D sample, i.e. D° 
and D° (-). 

The mixing parameters x" and y", 

( x" \ ( coa6 Kmr sin 6 Kmr 
\y" ) V - sm ^7T7r cos 5^ 

may be measured by means of a time- dependent Dalitz plot analysis of the 
three-body WS events D° — > K + 7t~tt°. The strong phase that rotates the 

mixing parameters x and y is defined here as = arg ( ^{b^k^-] 

5. -0(3770) — > D°D° reaction: exploiting the quantum correlations of the D°D° 
production at threshold, it is possible to constrain the physics parameters by 
measuring x 2 , y, Rd, 2^R d cos 5, 2y/R D cos 5. 
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4 HFAG results for D°-D° mixing and CPV 



HFAG uses a global fit to determine world averages of mixing parameters, CPV pa- 
rameters, and strong phases. The fit uses 38 observables taken from measurements 
of D° -> K+fu, D° ->■ K + K~, D° -> vr+vr", L>° ->■ iT+vr", D° -> K+tt-tt , 
D° — )■ i^57r + 7r _ , D° — > K$K + K~ decays and from double-tagged branching fractions 
measured at the ?/>(3770) resonance. Correlations among observables are accounted 
for by using covariance matrices provided by the experimental collaborations. Errors 
are assumed to be Gaussian, and systematic errors among different experiments are 
assumed uncorrelated unless specific correlations have been identified. An indepen- 
dent log-likelihood fit which accounts for non-Gaussian errors has been used as a 
control check, and equivalent results have been obtained. 

Averages for the parameters Rm, Vcp an d A r are calculated and then provided 
as input to the global fit. The observable Rm is calculated from D° — > K + £~v 
decays [201 12H 1221 [23] , and the average value is R M = (0.013 ± 0.027)%. The inputs 
used for this average are plotted in Fig. [TJ The observable yep is calculated from 
D° — > K + K~ /tt + tt~ and from D° — > KgK + K~ decays while is calculated from 
D° — > K + K~ /tt + tt~ decays only. The inputs used for these averages are plotted in 
Fig. H The average values are y C p = (1.064 ± 0.209)% and A T = (0.026 ± 0.231)%. 

The D°^K + n- observables used are from Belle [M], BABAR [1], and CDF ©; 
earlier measurements have much less precision and are not used. The observables from 
D° K° s 7T+7T- decays for no-CiV are from Belle [23] and BABAR [26] , but for the CPV- 
allowed case only Belle measurements [25] are available. The D° — > K + tt~tc° results 
are from BABAR [3], and the ^(3770) ->■ D°D° results are from CLEOc [27J. Time- 
integrated measurements of CP-violating asymmetries {Acp) in the D° — > K + K~ 
and D° — > -k + ti~ modes are from BABAR |29j and Belle [30J. The difference of the 
CP- violating asymmetries defined as AAqp = Aqp(K + K~) — Acp{^ + ^) are from 
LHCb [2] and CDF [3T]. These measurements reported evidence for direct CPV with 
a statistical significance of 3.5<r and 2.7<r, respectively. 

Three types of fit are produced by HFAG with different CPV assumptions: 

1. No CPV: in this fit it is assumed that CP is conserved and the CPV parameters 
Ap,, Ak, An, \q/p\-l and (ft are fixed to zero. 

2. No direct CPV: the CPV parameters Ad, Ak, A n are fixed to zero. In this 
case the relation [32], [33] tan0 = (1 — \q/p\ 2 )/(l + \q/p\ 2 ) x (x/y) is satisfied, 
and this reduces the four independent parameters (x, y, \q/p\, 4>) to three. The 
independent parameters used in the fit are x 12 = 2\M l2 \/T, y 12 = T 12 /T, and 
<p 12 = Arg(M 12 /r i2 ), where M 12 and T 12 are the off-diagonal elements of the 
D° — D° mass and decay matrices, respectively. 
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HFAG -charm 



E791 1996 



CLEO 2005 



BaBar 2007 



Belle 2008 



World average 



H 



H 



0.110 + „ 



0.160 ± 0.290 ± 0.290 ' 



0.004 + „ 



0.013 ± 0.022 ± 0.020 ' 



0.013 ± 0.027 % 





-0.2-0.1 0.1 0.2 0.3 0.4 0.5 



R M (%) 



Figure 1: World average value of Rm from Ref. [19], as calculated from D° — > K + £ v 
measurements [20112711221125]. 



3. CPV -allowed: where all the parameters for mixing and CPV are floated in the 
fit. 

All fit results are listed in Table HJ The total \ 2 is 35.6 for 37— 10 = 27 degrees of 
freedom; this corresponds to a confidence level of 12.4%. The resulting la-5cr contours 
are shown in Fig. [3] for the CP-conserving case, in Fig. H] for the no-direct-CPK case, 
and in Fig. [5] for the CPK-allowed case. For the CPK-allowed fit, the no- mixing point 
(x, y) = (0,0) is excluded at a confidence level of 1.28 x 10 -24 corresponding to a 
statistical significance of 10.2cr. The parameter x differs from zero by 2.7cr, and y 
differs from zero by 6.0a. In the (\q/p\,4>) plot, the point (1,0) is within the la 
contour; thus the data are consistent with CP conservation in mixing and in the 
interference between mixing and decay. 
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World average 



0.732 ± 2.890 ± 1.030 % 
3.420 ± 1.390 ± 0.740 % Belle 2007 

-1.200 ± 2.500 ± 1.400 % 
-0.500 + 1.000 ± 0.800 % BaBar 2007 
1.310 ± 0.320 ± 0.250 % 
0.110 ± 0.610 ± 0.520 % LHCb 2012 
1.160 ± 0.220 ± 0.180 % 
0.550 ± 0.630 ± 0.410 % 



1.064 ± 0.209 % 



-4 -3 -2 -1 1 2 3 4 5 



World average 



0.010 ± 0.300 ± 0.150 ' 



0.260 ± 0.360 ± 0.080 % 



-0.590 ± 0.590 ± 0.210 ' 



0.026 ± 0.231 % 



y CP W 
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Figure 2: World average value of y CP and A-p from Ref. [19], as calculated from 
D°^K + K-/n+ir-. 




Figure 3: Two-dimensional contours for the mixing parameters {x,y), for no CFV . 
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Figure 4: Two-dimensional contours for theory parameters (x 12 ,y 12 ) (top left), 
(x 12 ,0 12 ) (top right), and (?/i 2 ,0i 2 ) (bottom), for no direct CPV. 
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Figure 5: Two-dimensional contours for the parameters (x,y) (left) and (\q/p\,<j)) 
(right), allowing for CFV. 



Table 1: Results of the global fit for different assumptions concerning CFV. 



Parameter 


No CFV 


No direct CFV 


CFV-allowed 


CFV-allowed 95% C.L. 


x (%) 


n 65 + 018 

U - D0 -0.19 


0.62 ±0.19 


63 + 019 

U.OO _ .20 


[0.24, 0.99] 


y{%) 


0.73 ±0.12 


0.75 ±0.12 


0.75 ±0.12 


[0.51, 0.98] 


5(°) 


Z1 - u -11.0 


99 9 +9.9 
-H.2 


99 i +9.7 
z ^- 1 -11.1 


[-2.6, 40.6] 


R D (%) 


0.3307 ±0.0080 


0.3305 ± 0.0080 


0.3311 ±0.0081 


[0.315, 0.347] 


A D (%) 


0.0 


0.0 


-1.7 ±2.4 


[-6.4, 3.0] 


\q/p\ 


1.0 


1 04 + 07 

1-U1 _ .06 


0.88183 


[0.59, 1.26] 


<M°) 


0.0 


-2.02t|;?I 


-io.i iL 5 


[-27.4, 8.7] 


W (°) 


-i 7 o +21.7 
i '- 8 -22.8 


19.4™ 


i q q +21.8 

J-f.O _22.9 


[-26.3, 61.8] 


A„ 


0.0 


0.0 


0.36 ±0.25 


[-0.13, 0.86] 


A K 


0.0 


0.0 


-0.31 ±0.24 


[-0.78, 0.15] 


x 12 (%) 




0.62 ±0.19 






Vl2 (%) 




0.75 ±0.12 






0ia(°) 




-6.5 







10 



5 New results from the B factories 



5.1 Measurement of D°-D° mixing and search for indirect 
CP violation in D° — »■ K + K~ and D° — »■ 7r + 7r~ decays 

The BABAR [34] and the Belle experiments [35] have recently presented updated results 
for the measurements of the mixing parameter yep and the CP violation parameters 
AY for BABAR and Ar for Belle. The definitions of AY and Ar are the following: 

"P+ I 1- -j- 

AF = ~ , A r = Lf^r and AF = (1 + y CP )A r , (17) 
21 r + + t+ 

where r + = l/r + (r + = 1/T + ) are the effective lifetimes for D° (D°) decaying to 
the CP-even final states K + K~ and 7r + 7r~. In principle the parameters yep and 
AY depend on the final state /, as indicated in Eq. ( ITS]) and Eq. f[T9~j) below, when 
accounting for direct CP violation. The BABAR and Belle analyses assume CP con- 
servation in the decay (i.e. A n = 0, Ak = 0), and neglect terms of O(10 -4 ) in the 
expressions below, that are beyond the present experimental sensitivity: 

f 1/ 1 

2/£p = ycos(f> f + -{A M + A f )xsm(j) f - -A M A f y cos(p f (18) 

l/cp -2/ cos + sin 

1 1 

AY } = —xsm<fif + ~(A M + Af)ycos<f)f + -A M AfXsin<fif (19) 

^ AF ~ — xsin0 + -A M y cos <p. 

Hence, the parameters no longer depend on the final state /, and so can be averaged 
over the K + K~ and 7r + 7r~ modes. 

The measurements are based on the ratio of lifetimes extracted simultaneously 
from a sample of D° mesons produced through the flavor-tagged process D* + — > 
D°7T + , where the D° decays to K~tc + , K~K + , or 7r~7r + ; BABAR uses the additional 
samples of untagged decays D° — > K~n + and D° — > K~K + for the measurement 
of yep- These have about 4 times the statistics of the corresponding flavor-tagged 
sample, but have lower purity. 

The main selection criteria require that the center-of-mass momentum of the D° 
(p*) be greater that 2.5 GeV/cH, that D° daughter tracks be identified as kaons 
and pions, and require for the flavor-tagged sample that the reconstructed Am = 



•i-Note that this definition for Ay has different sign convention with respect to BABAR results 
published previously [21 0] . 

§Belle requires p* > 3.1 GeV/c for the data collected at the T(5S') mass peak. 
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Avr< P (%) A(AY) (%) 


Fit region 
Signal model 
Charm bkg 
Combinatorial bkg 
Selection 


0.057 0.022 
0.022 0.000 
0.045 0.001 
0.079 0.002 
0.059 0.054 


Total 


0.124 0.058 



Table 2: Systematic uncertainties for yep and AY at BABAR. The total is the sum- 
in-quadrature of the entries in each column. 



m(D* + ) — m(D°) be close to the value of 0.1455 GeV. The selection criterion on p* 
is used to reject D mesons from B decays, and to improve the signal significance. 

The proper time t and proper-time error a t are obtained from the reconstruction 
of the 3-dimensional flight length (L) and the momentum of the D° (p) according to 
the relation t = m/\p\ 2 L -p, where m is the nominal mass of the D°. The flight length 
is reconstructed by means of a kinematic fit to the decay vertex and production vertex 
of the D°, the latter being constrained to originate within the e + e~ collision region. 
The typical transverse dimensions of the luminous region of the PEP-II collider are 
about 100 fim in the x direction and 7 /im in the y direction. The most probable at 
value is about 40% of the nominal D° lifetime, and only candidates with <j t < 0.5 ps 
are retained for the fit in the BABAR analysis. 

The lifetimes of the CP-even modes K~K + , 7t~tt + are compared to that of the 
CP-mixed mode K~ir + in order to measure yep, which is proportional to the ratio 
of the lifetimes, and AY (A?) which is proportional to the difference of the effec- 
tive lifetimes of D° and D° into CP-even modes. BABAR measures yep = [0.72 ± 
0.18(stat) ± 0.12(syst)]% and AY = [0.09 ± 0.26(stat) ± 0.06(syst)]% using a data 
sample corresponding to an integrated luminosity of 468 fb~ l [33]. Belle measures 
y C p = [1.11 ± 0.22(stat) ± 0.11(syst)]% and A T = [-0.03 ± 0.20(stat) ± 0.08(syst)]% 
using a data sample corresponding to 976 fb' 1 [35] . 

The systematic uncertainties on yep and AY are reported in Table [2] for BABAR, 
and on yep and Ar are reported in Table [3j for Belle. The total systematic uncer- 
tainties are comparable between the two experiments. It is worth noting that the 
Belle experiment quotes a systematic error due to the silicon vertex detector (SVD) 
misalignment that is is negligible in the case of the BABAR experiment. The measure- 
ment of the proper time average value in Belle shows a dependence on the cosine of 
the polar angle in the e + e~ center-of-mass (cos#*) that is not properly reproduced by 
Monte Carlo simulation. The discrepancy between data and Monte Carlo is different 
for the data collected with the two vertex detectors that have been used during the 
running of the experiment. In order to minimize the systematic error due to this 



12 



C!^ tPPTYTV 


Av^p (%) AAr (%) 


Acceptance 

SVD misalignments 

Mass window position 

Background 

Resolution function 


0.050 0.044 
0.060 0.041 
0.007 0.009 
0.059 0.050 
0.030 0.002 


Total 


0.11 0.08 



Table 3: Systematic uncertainties for yep and Ay at Belle. The total is the sum-in- 
quadrature of the entries in each column. 



effect, Belle has performed the measurement by dividing the sample into different 
intervals of cos#*. 

The yep measurements report evidence for D°-D° mixing with a significance of 
3.3cr in the case of BABAR (most precise measurement to date) and 4.5a in the case 
of Belle. The measurements of AY from BABAR, and Ay from Belle are compatible 
within error and both are consistent with no CFV . The new measurements supersede 
the previous BABAR (2JH] and Belle results [5]. The updated HFAG averages, including 
the new results, are y C p = (0.866 ± 0.155)% and A r = (-0.022 ± 0.161)%. The 
comparison with the previous HFAG average values yep = (1.064 ± 0.209)% and 
Ay = (0.026 ± 0.231)% indicates significant improvement in precision and a lower 
central value for y CP . 

6 Summary 

In summary, HFAG provides averages for the D°-D° mixing and CP-violating pa- 
rameters by combining results from the following experiments: BABAR, Belle, CDF, 
CLEO, CLEOc, E791, FOCUS, and LHCb. The average values from the CiV-allowed 
fit for the mixing parametes are x = (0.63l ; 2 o)%, V = (0.72 ± 0.12)%, and for the 
CP-violating parameters are \q/p\ = 0.88^11, = (— lO.llgJ) . These average val- 
ues do not include the new BABAR and Belle results presented here. The parameter x 
differs from zero by 2.7a, and y differs from zero by 6.0cr. The no-mixing hypothesis, 
(x = 0,y = 0), is excluded with a statistical significance of 10.2<r. The LHCb and 
CDF experiments have obtained first evidence for direct CFV in D° decays, while 
there is no evidence for CFV in mixing (\q/p\ ^ 1) and in the interference between 
mixing and decay (0^0). 

New results for D°-D° mixing and CFV in the lifetime ratio analysis of the 
transitions D° — > K + K~ /tt + tt~ to the transitions D° — > K~n + have been pre- 
sented. BABAR measures y CP = [0.72 ± 0.18(stat) ± 0.12(syst)]% and AY = [0.09 ± 
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0.26(stat) ± 0.06(syst)]% using a data sample corresponding to an integrated lu- 
minosity of 468 /ft -1 . Belle measures y CP = [1.11 ± 0.22(stat) ± 0.11(syst)]% and 
A r = [— 0.03±0.20(stat)±0.08(syst)]% using a data sample corresponding to 976 fb^ 1 . 
The updated HFAG averages, including the new results, are yep = (0.866 ± 0.155)% 
and A r = (-0.022 ± 0.161)%. 

HFAG averages for D°-D° mixing and CFV parameters are in agreement with 
the SM. Recent results for direct CFV in D° — > K + K~ /it + it~ decays report a larger 
asymmetry value than the SM expectation. However, given the present knowledge of 
the charm system, this cannot (yet) be considered a clear signal of physics beyond 
the SM. 
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